Fermat's principle in quantum gravitational optics 

N. Ahmadi a *, S. Khoeim-Moghaddam 6 '^ and M. Nouri-Zonoz a,c * 
a Department of Physics, University of Tehran, 
North Karegar Ave., Tehran 14395-547, Iran. 
b Department of Physics, Sharif University of Technology, 
P Box 11365-9161 Tehran, Iran. 
Institute for studies in theoretical physics and mathematics, 
P Box 19395-5531 Tehran, Iran. 

Abstract 

Interactions incorporating the vacuum polarization effects in curved backgrounds modify the null 
cone structure in such a way that the photon trajectories would not be the space-time geodesies 
anymore. The gravitational birefringence introduced as a direct consequence of these effects, will 
allow shifts in the photon velocities leading to polarization dependent superluminal propagation. 
Taking these effects into account we study Fermat's principle in the context of the 1+3 (threading) 
formulation of the space-time decomposition. We find an expression for the modified spacetime 
refractive index and show it is proportional to the light cone correction to the first order. Conse- 
quences of this modification on polarization sum rules and spatial light paths are considered. 
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I. INTRODUCTION 



Quantum Gravitational Optics (QGO) is the term coined for the studies incorporating 
photon propagation in curved space-times in the presence of QED vacuum polarization in 
the lowest (leading) order . The quantum corrections introduce null cone modifications 
which in general alter the propagation characteristics of photons leading to cases in which 
they could travel at speeds greater than unity. It is expected that the more complex the 
space-time is, the more complex should be its curvature coupling to these quantum effects. 



In line with the above expectation, this phenomenon 
Robert son- Walker and gravitational wave backgrounds 



ras been studied in Schwarzschild, 
l| as well as in Reissner-Nordstrom 



and Kerr black hole space-times 0, 4]. It is found that except for the case of de Sitter space 
there may be certain directions and polarizations in which photons travel at speeds greater 
than unity. It is also shown that in the case of static topological black holes, the velocity 
shift of photons is the same as in the Reissner-Nordstrom case and the vacuum polarized 
photons are not sensitive to the asymptotic behavior or topological structure of space-time 
while for radiating topological black hole backgrounds, the light cone condition and velocity 
shifts depend on the topological structure |5(. In a recent paper, taking the vacuum polar- 
ization into account the propagation of a bundle of rays in a background gravitational field 
has been studied through the perturbative deformation of the Raychaudhuri equation {(J. 
However, there are other important phenomena predicted in the framework of general rel- 
ativity (GR) while semiclassical corrections such as the QED interactions mentioned above 
are taken into account. Despite the fact that QGO corrections are expected to be small, the 
comparison between the classical (GR) and semiclassical (QGO) theories can be conducted 
on the basis of higher order effects which deserve attention from both the theoretical and 
observational points of view. In the realm of photon propagation, the observational testbeds 
are usually provided in the context of gravitational lensing through the impressive devel- 
opment of technical capabilities. Hence a special treatment of lensing phenomena could be 
provided through the consideration of these quantum effects. Following a particular version 
of Fermat's principle for stationary space-times, we obtain an equation for the spatial light 
paths which demonstrates why, for a QED photon, a curved background acts as if it is an 
anisotropic medium with its own polarization dependent refractive index. 

In what follows, first we give a brief account of QGO. Then, based on Fermat's principle 



2 



in the context of 1 + 3 (threading) formulation of space-time decomposition, it is shown how 
the basic modifications introduced by the theory lead to the corrections on the photon spatial 
trajectories and the spacetime refractive index. After introducing an intriguing expression 
for an effective local refractive index in section 3, it is shown in section 4 that its first order 
correction is proportional to the null cone modification thereby leading to two new theorems 
on a)-polarization summed shift of refractive indices and b)-photons travelling normal to 
the horizon. In section 5 using the results of Tsai and Erber ?J on the refractive indices 
of photons of arbitrary frequency propagating in a background magnetic field, gravitational 
and electromagnetic birefringence are compared. In section 6 the equation of spatial light 
paths in the presence of the above QGO effects is given. The results are summarized and 
discussed in the last section. 

II. VACUUM POLARIZED PHOTON PROPAGATION IN A CURVED BACK- 
GROUND 

Vacuum polarization is an effect in which the photon exists as a virtual e + e~ pair for a 
short time. This virtual transition assigns photons with an effective size of 0(\ c ), where 
A c is the Compton wavelength of the electron The photon propagation will be affected 
by gravity if the scale of space-time curvature L is comparable to A c . Ignoring the vacuum 
polarization, the equivalence principle leads to the photon propagation at the speed of light 
on the space-time null geodesies. But when the pair production is taken into account, the 
equivalence principle is violated in such a way that the photon trajectories are modified and 
superluminal (photon) propagation becomes a possibility. All this could happen without 
necessarily breaking down the causality Considering the vacuum polarization, the first 
order correction to the electromagnetic action 




(1) 



is given by 



W, 




x(-g)*(aRF^F>" + bR^F^F^ + cR^F^ F aT + dD^D v F v a 



) (2) 



in which 



a = 



1 a 
144T 



b 



13 a 
3G0n 



c = 



1 a 
360n 



d 



1 a 
30n 



(3) 
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are perturbative coefficients of O(a), m e is the electron mass and a is the fine structure 
constant. Expression (j2J) is called Drummond-Hathrell action [1]. The first three terms show 
the influence of curvature while the last one, represen ting the off-mass-shell effects (in the 
vacuum polarization), exists even in the flat spacetime 14|. The effective equation of motion 
can be obtained by varying W = W + W\ with respect to A v as follows, 

D (l F ,a/ - — [2bR^D"F av + AcR^ aT D»F aT ] = (4) 

in which only terms of first order in a and (A c /L) 2 are kept. Now we use the geometrical 
optics approximation, in which the electromagnetic field is written as a slowly-varying ampli- 
tude and a rapidly-varying phase through a small parameter e, i.e, (A fi + ieB li + - ■ ■ ) exp(z^), 
where d is a real scalar field. In geometric optics, the wave vector is defined as the gradient 
of the phase, i.e, k^ = d^d and the polarization vector is introduced through = Aa^, 
where A is the amplitude and a M is the polarization vector. Rewriting the equation of motion 
(J3J) as an equation for the polarization vector we have |8j, 

2b u 8c 

k a a H - iRnuik^k a a k^k^a )] H ■^Ra aVT V , k a =0, (5) 

From which the light cone condition at 0{e~ 1 ) is given by 

2b 8c 

k H -^Riivk^k -Ruaurk^k a a =0, (6) 

in which a M a M = — 1. Since this effective equation of motion is homogeneous and quadratic 
in k^ it can be written as 

g^ u k^k v = o, (7) 

in which Q^ v could be defined as an effective (physical) metric whose null geodesies represent 
the modified photon trajectories. In other words at this order of approximation there is no 
dispersion and consequently the phase and group velocities are equal and given by 

k° 

v ph = (8) 

From (j7]) it could be seen that there is a possibility for the modified null cone to lie outside the 
original light cone leading to superluminal photon velocities as well as modified trajectories. 
It is this feature of QGO that will be the main focus of our study through Fermat's principle 
in the forthcoming sections. 
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III. FERMAT'S PRINCIPLE IN QGO 



One way to look at the nature of the modified null trajectories in QGO is through 
the fundamental definition of null rays via Fermat's principle. Within the assumptions of 
geometrical optics, the trajectory of a light ray 7 from a source P to an observer on his 



world line C, in a space-time 
version of Fermat's principle 



yM,<7 M i,), can be characterized by the following relativistic 
9J: Signal's arrival time is stationary with respect to the first 
order variations of 7 within the set of smooth null curves from P to C . 
In stationary space-times Fermat's principle takes a version of particular interest. To invoke 
this version of the principle, we employ the 1 + 3 (threading) formulation of space-time 
decomposition in which the space-time metric can be written in the following form |lo| 

ds 2 = h (dx° - Adx*) 2 - all 2 (9) 

where 

A i = gi = , h = g 00 

9oo 



and 



dl 2 = ^ijdx l dx > = I —gij + J dx % dx 3 , i,j = 1, 2, 3 

goo 



is the spatial distance in terms of the 3-dimensional metric 7^ of the three dimensional 
space. Using the identity g^ v g va = 5% it could be shown that the following relations also 
hold 

>yU = -g<S , g l = -g° l , g™ = —-g t g* (10) 

#00 



Using the classical version of Fermat's principle (i.e at O (a )) [lOj, we obtain its QGO 
version for the propagation of null rays in an stationary gravitational field, taking into 
account the above vacuum polarization effects. Fermat's principle states 

5 J kidx 1 = 0, (11) 

where the integral is taken along the spatial projection of the physical rays which are to be 
varied with fixed endpoints. Rewriting the light cone equation (jBJ) in the form, 

= \ (%cR^ UT a T a a - 2bR ia/ ) , (12) 

7Tt„ 
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we obtain 

h (k° - gik*) 2 - -ufkV = A 00 k°k + 2A 0i ~k°k i + A^P. (13) 

In which k stands for the modified wave vector to the first order in a. Noting that the 
vector k % must have the direction of the vector dx 1 for the classical photon (zeroth-order 
of the above equation), i.e, k l = -^^ff we find its correction of O (a) by substituting the 
classical values for k l and k° in the right hand side of (fl3|) . and the conjectured modified 
wave vector k % = + B l in its left hand side so that it reduces to 



k ** + Bl){ Kd^_ + Bj) = AoQ{kr + 2Aoik oh>d^ + ^tf**^ (14) 



k_l_ 

h lijK Vh dl ' " /v v^ dl 
where we have used the fact that 



Vh dl 



h dl dl 



k = g 0i k l = h(k° - g t k l ) = h(k° - 9i V) 



(15) 



is the constant frequency along the ray. After a few lines of simple manipulations it is found 
that 

(16) 



B = ~2ko { Aoo{k ] H + 2A ^ 7T + ~h Akn ~dT 



and consequently substituting for k° from (fl5|) we end up with 



~. = k^dx^ f _ A q d2 



2A 0j ^ l D + A kj y> % 



a dx k 
~~dT 



(17) 



Vhdi y 2 - ) 2 Vh\ j 

in which D = ^ + <?i~|~- It is not difficult to double check that this vector satisfies equation 
( TTBl up to the first order in a. To obtain the expression for ki we use the following relation, 

k i = gH IM = -g i k -^k j (18) 

which is a direct consequence of the relations given in ( flOl) . Therefore 

h = -7»j (k J + k gA 

d >■■' ( 4,„, 



-ho 
-k 



9i + lij 
9% + lij 



1 dx^ 




Vh dl 




1 dx^ 




dl 





— 7ij (2A 0kl k W + A kl ^— 



2y/h 



dl 



D- 



1 ( . dx 
( 2A 0iJ D + A ki 



2Vh V 



(19) 



Multiplying by dx 1 and dropping the constant factor fc we en d up with the following modified 
version of the Fermat principle, 

D . , . , 1 . . A dx j 
~7h 
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J A(l - ^D 2 ) + (ft - -^4ot)<fc* - -^=^dx*=r- = 0. 



2y/h 



dl 



(20) 



In a static background (for which D = -^=) we have 

M 1 - -sr)-p - ~4 ^f= — = 0. 21 

J 2h y/h h 2Vhdl | p. 

Comparison of (|20|) with Fermat's principle in its classical form, i.e S J ndl = |9|, Il0| . 

indicates that the trajectory of a light ray in a gravitational field is determined in the same 

way as in an inhomogeneous refractive medium with n as its effective local index of refraction 

given by 



n = 7K {1 ~ ^° 2) + i9i ~ 7K Aoi)el ~ 27^' eV ' (22) 



n = E°- ^^A^E^E". (24) 



in which e* = ^ is the tangent vector to the ray. The equations f[2"Uj) . (12"T1) and ( J22D reduce 
to their corresponding classical counterparts at zero a limit. An alternative representation is 
given by introducing the four vector E^ = =^-, which is null with respect to the background 
metric and its components are (_D,e J ). The covariant components of the tangent vector e* 
form a three-dimensional vector in the space with metric 7^ and correspondingly the square 
of this vector is equal to one. 

e» = j^e? e 2 = 7y-eV = 1 (23) 

With such a definition, the effective local index of refraction is expressed as 

1 

2V7T 

This is obviously a result more general than the standard gravitational result discussed 
in 10j where only the effect of the gravitational field (in the context of 1 + 3 approach), 
as a dispersive medium, on photon propagation is considered. A comparison between the 
classical wave vector = (k°, -^e*j and E^ shows that = -^E^ i.e, these two vectors 
are parallel . Therefore we have 

5n = --^A^E" = -T^A^kr = -||(^ 2 ), (25) 

where 5n and 5k 2 are refractive index and null cone modifications respectively. The 
has been studied for different directions and polarizations in many curved backgrounds 
So evaluation of 5n reduces to the multiplication of 5 k 2 by the local factor — 



atter 

Ik 



IV. POLARIZATION SUM RULES AND HORIZON THEOREM 



In the previous section it is shown that Sn is proportional to 5k 2 , using ([25]) we can state 
two theorems parallel to polarization sum rules and horizon theorem in J8J. 
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A. Polarization sum rules 



For Ricci flat space-times, the sum over the two physical polarization of the refractive 
index is zero i.e, Sn = 0, but for non Ricci flat space-times satisfying the Einstein field 

pol 

equations, 

E 5n = ^^ b + 2c ) T ^ kr ( 26 ) 

pol 

Substituting for coefficients a and b ,we arrive at 



E^=i4^V r ' (27) 
^— ' 45 wr ki 

pol u 

where T^ u is the energy momentum tensor. Since in the weak field perturbative expansion, 
we can take to be null on the right hand side of equation (I27p . the null energy condition 
implies T^ v k^k u > 0. So assuming the null energy condition, we always have ^ £n > 0. 

pol 

This relation is consistent with the polarization summed velocity shift, Y^Sv < 

pol 

B. Horizon theorem 



ii 



In 



,8 



it is shown that for spacetimes with an event horizon, photons with momentum 
directed normal to the horizon have velocity equal to c, i.e. the light cone remains k 2 = 
independent of photon polarization. Using fl25l) yields that such (QGO) photons also have 
local refractive index equal to E° , independent of their polarization. 



V. ELECTROMAGNETIC ANALOGY 

With a similar discussion, one can find the modified refractive index corresponding to the 
modified light cone condition in an arbitrary anisotropic but homogenous electromagnetic 
field in flat space-time. The description of such a photon based on the one-loop action 
can be given by an expansion in powers of which is a dimensionless parameter defined 
as the ratio of the background scale to a typical mass in the theory. Starting from the 
Euler-Heisenberg [kJ effective action 



T = j dx [-^F^ + -L (z (F^f + yF, v F Xa F^F»°) 



(28) 



with z = — -^a 2 and y = ^a 2 , we derive the equation of motion 

Analogous to the equation (jSJ) in the gravitational case, using similar arguments the geo- 
metric optics approximation gives 

W" - — jF^Fxak^f^ - -4 {F» x F Xa k,r + F» x F™k,f Xa ) = 0, (29) 
m 4 m 4 v ' 

where F^ v is the background electromagnetic field strength. Rewriting the above equation 
as an equation for the physical polarization vector, a^a^ = — 1, we have 

k 2 + — A (4* + y) F^F Xa k»k x a»a° - ^F^F Xl/ k^k x = (30) 
m 4 m 4 p 

For a polarization state orthogonal to the plane spanned by the photon momentum and 
magnetic field direction, in which F^af 1 = 0, the light cone condition is 

k 2 = %F x F vX F^k v = -^T^k v (31) 
m 4 p m 4 

Here in the second equality with the classical energy-momentum tensor expression, T^ v = 
—F^ x F uX + \r]^ u F Xa F Xa , the higher order corrections are neglected. For the second state 
coplanar with the above plane, it is given by 

k 2 = -—(2z + y)T^k» (32) 

where in both cases (l3Tj) and (1321) . the null cone dependence on the energy momentum tensor 
is similar to that of the gravitational case and the polarization dependence appears through 
the coefficient of T fJiU k tJ, k 1 ' (or what we have called A^ v in equation (|12p ). Setting h and g a 
equal to one and zero respectively in defination of components and using equation fl24l) . 
the local refractive index for each polarization is given by 

1 4y 

Kq m 4 

n 2 = 1 + i^r (?z + y) T^k v (33) 
Kq m 4 



Again notice that if T^ v k^k u > 0, both polarizations have refractive indices larger than one. 
It is also notable that the above expressions agree with the results of j3] in the weak field 
and low frequency limit. Now taking the polarization sum and substituting for the explicit 
values of coefficients y and z, we find 

pol u 



This should be compared with eq. (127)) for the gravitational case. Notice that for electromag- 
netism, the sum over the one- loop correction is proportional to a 2 , essentially because the 
lowest order terms in the Euler-Heisenberg action are of 4th order in the background fields, 
i.e. 0(F A ), compared to the 3rd order terms (i.e. O(RFF)) in the Drummond-Hathrell 
action 



VI. SPATIAL LIGHT PATHS 



The stationarity of the photon arrival time under the first order variation of 7 spatially 
translates into the stationarity of its arbitrarily parametrized spatial arc length 



dl 



5 I n(x i ,e i )—d\ = 5 \ n(x i ,&)y/¥E i d\ = Q, 



(35) 



in which A is a general parameter and e l 



c co- 
variation reduces to the spatial light path equation 

d „ 



For an stationary spacetime the above 



d ( dn 
dx iU ~ dX yd^ 



(36) 



in which n = n(x l , e l ) yWei. Now choosing A for the extremal path to be proportional to the 
spatial arc length I, i.e A = al, we end up with the following equation for the light paths 



^ n ^ )= dJ 



d { d(ny/e*el] 



(37) 



Now for the classical part of n, i.e in the absence of the quantum corrections, we have 



d_ f d{E ^eej 
Jl [ de* 



d_ 
dl 



-jjZi + 9i+ {9je J ) e, 



(38) 



After some manipulation we get 



d 



E u —e 
dl 



(Vjjj -j= + e 3 (digj - d j g l 



e x V x g 



(39) 
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where Vj_ = V— e .Vj denotes the projection of V onto the plane normal to the 
direction of the ray given by the unit tangent vector e . This equation describes how a 
classical light ray in a curved stationary background is deviated from a straight line in 
Minkowski space. Deviation is devided into two intrinsically different parts. The first 
term corresponds to an attraction towards the deflection center, the so called gravitoelectric 
force, whereas the second term is due to the gravitomagnetic force originated from the 
gravitomagnetic potential 9 11]- 



In the presence of the quantum interactions, the equation for the spatial light paths, eq. (!37I) . 
is given by 



0, E [) - 0, I ^LE'eA = | 
2y/h / dl 



d ^_ / ^jiv e^E" 



(40) 



where = A /lu y/e i ei and the correction term in right hand side i 



is 



d_ 

dl 



2Vh Vh \dei 



+ e . V 



V ( d 



Gathering all together, we have 



1 d 
n dl 



e . V 



V 

.4 



) E u} + 



A, 



A 



= ei [ e . V 
d ( d 



2Vh J 2Vh dl 



dl \de l J V de l J dl 



(41) 



^E"E V + 
Vh 2Vh J 

>'" ' —E { ^ | E u} + A>J,U 



e x V x g 



- ( #^ ET> + ( ±E*A ±E»> 
dl\de l J V de l J dl 



,(42) 



in which A^ L B u ^ = | (A^B U + A u B fl ) and the derivatives are given by the following relations 
in the same coordinate frame 







di [ -j= ) + e 3 digj, 



= M 

d /- -+ ( 1 \ - - , , i d , d 



(43) 



dl \ \JhJ J ' J dl ' dl 

The equation (1421) is a complicated combination of curvature and metric components which is 
formulated through the stationarity of the arrival time under the variations of the null curves 
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connecting the source to the timelike observer in the presence of quantum corrections. From 
the above calculations, compared to the classical path, the time delays are expected to be 
of order a . Therefore the natural question one should address is whether these corrections 
pile up as the travel distance increases. If the answer is yes, then a deviation much larger 
than the scale of the theory could be achieved. A more detailed discussion of this topic will 
be given elsewhere 12]. 



VII. DISCUSSION 



It is natural to expect that the inclusion of QED vaccum polarization effect on photons 
propagating in a background gravitational field should modify their trajectories in that field. 
This is actually part of the ongoing studies under the general title of quantum gravitational 
optics (QGO). To look at the roots of this trajectory modification we have studied Fermat's 
principle in QGO in the context of threading (1 + 3) formulation of spacetime decomposition 
which is best adapted to the employed version of the principle. Taking the above QGO ef- 
fects into account, it is shown there is an effective Fermat's principle for the propagation of 
a ray in a general stationary background through the introduction of an explicit expression 
for the effective index of refraction. Obviously not only QGO photons do not propagate 
along the shortest line in space but also the effective refractive index is itself a complicated 
combination of curvature and metric components as well as the quantum corrections con- 
sidered in the context of QGO. 

The modification of the refractive index is shown to be proportional to the null cone modifi- 
cation in the presence of QGO effects i.e 5k 2 . This fact is used to show that the polarization 
summed shift in the refractive indices is greater than or equal to zero. It is also shown that 
for photons with momentum directed normal to the horizon the refractive index being equal 
to the classical case is polaraization independent. 

Polarization dependent refractive indices for the electromagnetic analogue of the same prob- 
lem are calculated explicitly and compared with their gravitational counterparts in the weak 
field and low frequency limit. Finally, through the effective refractive index, the equation of 
spatial light paths in the presence of the above quantum corrections is found. 
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